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ABSTRACT 



Aims. To estimate distances to dense molecular cores that harbour Very Low Luminosity Objects (VeLLO) detected by Spitzer Space 
Telescope and to confirm their VeLLO nature. 

Methods. The cloud distances are estimated using near-IR photometric method. We use a technique that provides spectral classifica- 
tion of stars lying towards the fields containing the clouds into main sequence and giants. In this technique, the observed (J - H) and 
{H - Ks) colours are dereddened simultaneously using trial values of Ay and a normal interstellar extinction law. The best fit of the 
dereddened colours to the intrinsic colours giving a minimum value of then yields the corresponding spectral type and Ay for the 
star The main sequence stars, thus classified, are then utilized in an versus distance plot to bracket the cloud distances. The typical 
error in the estimation of distances to the clouds are found to be ~ 18%. 

Results. We estimated distances to seven cloud cores, IRAM04191, L1521F, BHRUl, L328, L673-7, L1014, and LI 148 using 
the above method. These clouds contain VeLLO candidates. The estimated distances to the cores are found to be 127 ± 25 pc 
(IRAM04191), 136 ± 36 pc (L1521F), 355 ± 65 pc (BHRlll), 217 ± 30 pc (L328), 240 ± 45 pc (L673-7), 258 ± 50 pc (L1014), 
and 301 ± 55 pc (LI 148). We re-evaluated the internal luminosities of the VeLLOs discovered in these seven clouds using the dis- 
tances estimated from this work. Except for L1014-IRS (L;„, = 0.15 Lq), all other VeLLO candidates are found to be consistent with 
the definition of a VeLLO (L,„, < 0. IL©). In addition to the cores that harbour VeLLO candidates, we also obtained distances to the 
clouds L323, L675, L676, CB 188, LI 122, LI 152, LI 155, LI 157 and LI 158 which are located in the directions of the above seven 
cores. Towards L1521F and LI 148 we found evidence of the presence of multiple dust layers. 

Key words. ISM: dust extinction - ISM: clouds - Stars: distances - Stars:low mass, brown dwarfs - infrared: ISM 



1. Introduction 

Some of the important physical parameters that are required to 
elucidate the still mysterious processes involved in the forma- 
tion and evolution of both molecular clouds and pre-main se- 
quence (PMS) stars depend crucially on the accurate determina- 
tion of their distances. But distances to most of the star form- 
ing and starless molecular clouds, especi ally those that are rel - 
atively isolated, are highly uncertain ( Hilton & LahuUaL Il995h . 
Furthermore, the need to determine accurate distances to cloud 
cores can be understood from the recent discoveries of embed- 
ded sources having luminosities less than Q.ILq called, very 
low luminous objects (VeLLOs), in a nu mber of dense cor es 
that wer e previou s ly con s idered as starle ss (lYoung et al^^ 2004. 



Kau ffmann et a 



2006c iLee et al. 



. .20051 iDunhamet ali 2006: Bourk e et al „ 
2009t iDunham et all l20i 0a). But, because the 



luminosity depends on distance squared, small errors in distance 
are effectively doubled. Therefore the above conclusions on the 
status of these objects as VeLLOs depend crucially on the dis- 
tances to the parent clouds, which are highly uncertain (in some 
cases by over a factor of 2). For example, the authors who car- 
ried out studies on the VeLLO in L101 4 have assum ed a distance 
of 200 pc to the cloud . But recentlv. iMorita et all (2006) have 
quoted a distance of 400-900 pc for L1014. Then the VeLLO, 



with luminosity of L ~ O.O9L0 (calculated using the assumed 
distance of 200 pc to the cloud), becomes 0.36-1 .8L0 and ceases 
to be a VeLLO. 

Dense cores were previously considered to be starless if 
they do not contain an Infrared Astronomical Satellite (IRAS) 
point source to a sensitivity of L ~ 0. 1 Lgid/ 140)^. But 
Spitzer observations of L1014, considered as a starless core, 
as a part of Spitzer Legacy progra m "From Mo lecular Cores 
to Planet Forming Disks" (or c2d; lEvans et all 12003 ) came 
as a surprise! LI 014 was found to contain an embedded 
source, L1014-IRS, with a very low l uminosity of L ~ 
0.09Le dYoung et all l2004t iBourke et all 2005; Huar d et all 
2006). Currently, the VeLLOs in fourteen more cores are 
identified in the full c2d sample jDunham et al.L |2008|) . Of 
these, five of them (including L1014-IRS) h ave been stud- 
ied in detail so fai-: IRAM 041914-1522 ('Dunham et all 
2006), L1521F-IRS (Bourke et al., 2006; Terebev et al., 200% 
L328- IRS (iLee et all I2009I). L673-7-I RS ( Dunham etaU 
2010ahandL1014 -IRS (lYoung et allf2004tlBourke et aU,2006[ 
Huard et alll2006h . Sources with such low lu minosities are hard 
to explain on the basis of the standard model of' Shuet al.l ( ll987h 
as the observed luminosities are found to be more than an or- 
der of magnitude lower than what is expected from the model 
assuming a spherical mass accretion onto a protostellar ob- 
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ject l ocated on the stellar/substellar boundary (iDunham et all 
l20Q6h . Thus, the discoveries of VeLLOs in dense cores can 
pose a formidable challenge to our understanding of star forma- 
tion. Recentlv. iDunham^ aO (|2010b) showed that the luminos- 
ity problem can be resolved and bring the model predictions in 
better agreement with the observatio ns if the episodic mass ac- 
cretion based on the simulation by Vorobvov & Basu (2006) is 
included in the evolutionary models of Young & Evans (2005). 
Then the sources with lowest luminosities are those that are ob- 
served in quiescent accretion states. 

A number of methods have been applied by various au- 
thors in the past to determine distances to clouds that are 
either associ ated with large comp l exes o r are relative l y iso- 
lated (e. g.. iGreenstein & Shaplevl 119371: iMcCuskevl 119391: 

Racinel flQ^sf InnttlipK Xj TTncnJ nOfiOl lRr>V Kr Mr-rnrtViJ 



Table 1. The J2000 coordinates of the cores as obtained from 
SIMBAD. 



lOottlieb & Upson', '1969": iRok & McCarthv' 
'1979; Tomitaetal., 1979: Snell, 1981; 
Reipurth & Gee, 1986; Cernis, 1990; Cernis & Straizvs, 1992; 
StraizysetaL .19921: Hilton & Lahulla. 1993; Kun&Prusti, 



1993MKenvon et allll9 94: Bourke et al., 1995a b; Corradi et al 
19971 iK nude & H0g, 1998; Peterson & Clemens, 1998; 
Nielsen et al., 2000; Franco, 2 002; Maheswar et al., 2004 
Loinard etal 



2007; 



2007 : iMaheswar et al. 



Alves & Franco, , _ „ - - - - 

2010^ ,Dzib et aL .2010,; .Knudel 120101: iKnude & Kaltchevai. 
2010h . Among the methods, the trigonometric parallaxes mea- 



sured using multi-epoch very long baseline array observations 
of radio-emitting young stars that are associated with the clouds 
give distance s to the corres p ondin g clouds with a prec i sion o f 
about 1-4% (iLoinard et all l2007l 120081: iTorres et all l2009h . 
However, this method can be applied only to those clouds 
that harbour radio-emitting young stars. Most of the other 
methods, mentioned above, are in general e xtremely tedious 
and d emanding in terms of telescope time. IMaheswar et al.l 
(12004 in an earlier study utilized optical and Two Micron all 
Sky Survey (2MASS) data to determine distance to Cometary 
Globule, CG12. One can combine the optical data from the 
Naval Observatory Merged Astrometric Dataset (NOMAD) 
and th e 2MASS data in the method given by Ma heswar et al.l 
(l200l to estimate distances to molecular clouds. The NOMAD 
database is basically a merge of data from the Hipparcos 
(Ferryman & ESA, 1997), Tycho-2 (H0g et al., 2000 ), UCA C-2 
tZacharias et al., .2004) and USNO-Bl (.Monet et al.r i2003') cat- 
alogues, supplemented by photometric data from the 2MASS. 
But the data available in this catalog is heterogeneous as the 
observations were carried at different epochs and also the 
accuracy of the photometry required by the method is high 
(< 0.1 mag) in BVR. Therefore it would be highly useful if 
we could develop a method that utilizes the vast homogeneous 
JHKs photometric data produced by the 2MASS, already 
available for the entire sky, to determine distances of the clouds. 
The n ear-IR photometr i c meth od presented in a previous pilot 
study (IMaheswar et all 1201 0') to estimate the cloud distances 
was such an endeavour. This method gives distances to clouds 
that are relatively closer (< 500 pc) with a precision of ~ 20%. 

In this paper we estimate distances to seven cloud cores 
namely, lRAM0419lQ, L1521F, BHR 111, L328, L673-7, 
L1014, and LI 148, which are identified to harbour VeLLO 
candidates dDunham et all i2008.) . using near-lR photometric 



' IRAM 04191-1-1522 is the name of the protostar found associated 
with a cloud core not listed in the catalogue by Lvnds (1962). This 
core was detected in the H'^CO^ J=l-0 survey by Onishi et al. (2002). 
The two cloud condensations identified by them towards this region are 
named as [OMK2002] 18a and [OMK2002] 18b. In this paper, here- 
after, we call the core that harbour IRAM 04191 + 1522 as IRAM04I9I. 



Core Identification 


Right Ascension 


Declination 




(° "' ■') 


r ' ") 


IRAM04I9I 


04 21 56.9 


+ 15 29 47 


L152IF 


04 30 50.3 


+23 00 09 


BHR 1 1 1 


15 42 20.0 


-52 49 06 


L328 


18 17 00.0 


-18 01 54 


L673 


19 20 54.4 


+ 11 13 12 


L1014 


21 24 05.9 


+49 59 07 


L1I48 


20 41 10.8 


+67 20 35 



method (IMaheswar et all l2010l) . Of the seven cloud cores, two 
of them, IRAM04191 and LI 52 IF are assocated with the Taurus 
molecular cloud complex for which jvery accurate distance mea- 
surements are available (ILoinard et alll2007l:ITorres et alll2007l 
12009 ). We included them in order to ascertain the reliability of 
the method. For the remaining five cores, most authors have used 
distances that are either guessed or assumed. The distance esti- 
mates available in the literature for each cloud core, prior to this 
study, are discussed along with our results in section 13.1! This 
paper is organized in the following manner: a brief description 
of the method and the data used to estimate cloud distances is 
given in section|2] The estimated distances of the clouds are pre- 
sented in section im The consequences of the new distance es- 
timates from this work on the status of the VeLLO candidates 
are discussed in section U!2] In section|4l we conclude the paper 
with a summary of our results. 



2. The data and the method 

Below we present a brief discussion on the methocfl We ex- 
tracted J, H and magnitude s of stars from th e 2MASS all- 
Sky Catalog of Point Sources (ICutri et all l2003h that satisfied 
the following criteria, (a) photometric uncertaint}0 cr < 0.035 
in all the three filters and (b) photometric quality flag of "AAA" 
in all the three filters, i.e., signal-to-noise ratio (SNR) > 10. We 
then selected the stars with their (J - Ks)< 0.75 to eliminate M- 
type stars from the analysis as unreddened M-type stars located 
across the reddening vectors of A0-K7 dwarfs make it difficult 
to differentiate the reddened A0-K7 dwarfs from the unreddened 
M-type stars. Also, the classical T Tauri stars (CTTSs), often 
found to be associated with the molecular clouds, occupy a well 
defined locus in the near infrared colour-colour (NIR-CC) dia- 
gram (Meyer et al. 1997) as shown in Fig. [T] which intercepts 
with the main sequence loci at (J - Ks)~ 0.75. Thus the criterion 
of {J-Ks)< 0.75 would allow us to eliminate most of the CTTSs 
as well. 

A set of dereddened colours for each star were produced 
from their observed colours by using a range of trial values of 
A V (0-10 mag) and the Rieke & Lebofsky (1985) reddening law. 
The choice of photometric uncertainty, cr < 0.035, ensures that 
the stars are selected from relatively low extinction regions of the 
cloud wher e the Rieke & Lebofsky (1985) redd ening law is most 
apphcable (IStraizvs et a"nil982HKandori et all l2003). The com- 
puted sets of dereddened colours of a star were then compared 



^ A more rigorous discussion on the errors and limitations of the 
method are presented in Maheswar et al. (2010) 

^ The photometric errors considered in this work in the 7, H, & 
magnitudes from the 2MASS database include the corrected band pho- 
tometric uncertainty, nightly photometric zero point uncertainty, and 
flat-fielding residual errors (Cutri et al. 2003). 
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Fig. 1. The (/ - H) vs. (H - K^) CC diagram drawn for stars 
(with Al/ > 1) from the region Fl towards IRAM04191 (see Fig. 
nil to illustrate the method. The solid curve represents locations 
of unreddened main sequence stars. The reddening vector for an 
AOV type star drawn parallel to the Rieke & Lebofskv (1985) 
interstellar reddening vector is shown by the dashed line. The 
locations of the main sequence stars of different spectral types 
are marked with square symbols. The region to the right of the 
reddening vector is known as the NIR excess region and corre- 
sponds to the location of PMS sources. Th e dash-dot-d ash line 
represents the loci of unreddened CTTSs (iMever et al.l Il997b . 
The line (J - Ks)< 0.75 is the upper limit set to eliminate M- 
type stars from the analysis as unreddened M-type stars located 
across the reddening vectors of A0-K7 dwarfs make it difficult 
to differentiate the reddened A0-K7 dwarfs from the unreddened 
M-type stars. The circles represent the observed colours and the 
arrows are drawn from the observed to the final colours obtained 
by the method for each star 



with the intrinsic colours of the normal main sequence stars. The 
intrinsic colours of the main sequence stars in t he spectral range 
A0-K7 were taken from lMaheswar et"al] (l2010l) . The best match 
giving a minimum value ofx^ then yielded the spectral type and 
Ay corresponding to that intrinsic colour. Once the spectral types 
and Ay values of the stars were known, their distances were es- 
timated using the distance equation, d (pc) - io*^'"'*^'-'+^"'**-)/^. 
Only those stars that were classified as dwarfs were considered 
for the determination of distances as the absolute magnitudes of 
giants are highly uncertain. The whole procedure is illustrated in 
FiglUwhere we plot the NIR-CC diagram for stars (with Ay > 1) 
chosen from the region Fl towards the direction of IRAM04191 
(see Fig.|4]i. The arrows are drawn from the observed data points 
to the corresponding dereddened colours estimated using the 
method. The maximum extinction values that can be measured 
using the method are those for AOV type stars (a; 4 magnitude). 
The extinction traced by stars will fall as we move towards more 
late type stars. 

In order to estimate distance to a cloud (not just by eye esti- 
mation), we first grouped the stars, classified as dwarfs, into dis- 
tance bins of bin width - 0.18 x distance. The centers of each 



bin were kept at a separation of half of the bin width. Because 
there exist very few stars at smaller distances, the mean value 
of the distances and the Ay of the stars in each bin were cal- 
culated by taking 1000 pc as the initial point and proceeded to- 
wards smaller distances. The mean distance of the stars in the 
bin at which a significant drop in the mean of the extinction oc- 
curred was taken as the distance to the cloud and the average of 
the uncertainty in the distances of the stars in that bin was taken 
as the final uncertainty of the determined distance of the cloud 
(see Figs.|2]and[3]). The vertical dashed line in the A^ vs. d plots, 
used to mark the cloud distance, is drawn at a distance deduced 
from the above procedure. The error in the mean values of Av 
were calculated using the expression, standard deviation /y/N , 
where is the number of stars in each bin. The Ay values were 
found to be uncertain by ~ 0.6 magnitude. The typical errors in 
our distance estimates for the clouds were found to be ~ 18% 
(iMaheswar et al.Ll2010t) . 

One of the prominent features of the molecular cloud struc- 
ture is that it is hierarchical. The extinction maps produced to- 
wards a number of regions in our Galaxy, using optical and 
near infrared techniques, have shown that the high-density fea- 
tures are inva riably contained insid e an extended low-densi ty en- 
velope (e.g., iDobashi et all 120051: iLombardi etalL l2010t) . Our 
method utilizes the background stars shining through the outer 
low-extinction (Ay < 4 magnitude) regions of the cores to esti- 
mate their distances. In case of the presence of a single dust layer 
in the line of sight, we expect a sharp increase in the extinction 
for the stars that are located behind that dust layer in an Ay vs. 
d plot. The distance at which the sharp increase in the extinc- 
tion occurs is considered as the distance to that dust layer. But if 
there exist more than one dust layer along the line of sight, we 
expect steep but step-like increase in the extinction at distances 
at which the dust layers are located. This, however, happens only 
if the foreground dust layers contribute a constant extinction val- 
ues. In case of patchy foreground dust layer(s), this method can 
give distance only to the first layer as the sharp step-like fea- 
tures get smeared by the non-uniform extinction caused by the 
foreground dust layei\ 

In a recent work, lKnudel(l2010h . using the 2MASS data, has 
presented a similar method to estimate distances to molecular 
clouds. Instead of using the main sequ ence stars in t he spec- 
tral range from A0-K7 as in our work, iKnudd (1201 Oh utilized 
three spectral groups: 0-G6 as primary sources, M4-T as sec- 
ondary sources and G6-M0 as teritary sources in the estimation 
of the cloud distances. The absolute magnitude in J-band were 
calibrated using Hipparcos stars. Following a rigorous statisti- 
cal approach, the jump in the extinction in the Ay vs. distance 
plot was obtained to determine the distance. In the pilot study, 
Knude (2010.) estimated distances to a number of nearby clouds 
that includes Chamaele on I and Lupus 3 which were included in 
IMaheswar et al.l (l2010l) also. The distance estimates in both the 
works are found to be in good agreeme nt with i n the error (the 
typical error in the distance estimates of iKnudd (1201 Ol) is found 
to be ~ 8%). 

3. Results and Discussion 

3.1. Distances to the cloud cores IRAM04191, L1521F, 
BHR111, L328, L673-7, L1014, and L1148 

We applied the method on seven dense cores, namely, 
IRAM04191, L1521F, BHRlll, L328, L673-7, L1014, and 
LI 148. We divided the fields containing the cores into small 
sub-fields to avoid complications created by the erroneous clas- 
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Fig. 2. The mean values of Ay vs. the mean values of distance 
plot for IRAM04191 and LI 52 IF produced using the the pro- 
cedure discussed in the ^to determine distances to the clouds. 
The error bars on the mean Ay values were calculated using the 
expression, standard deviation /y/N, where is the number of 
stars in each bin. 



sifications of giants into dwarfs. While the rise in the extinc- 
tion due to the presence of a cloud should occur almost at the 
same distance in all the fields, if the whole cloud is located at 
the same distance, the wrongly classified stars in the sub-fields 
would show high extinc t ion no t at the same but at random dis- 
tances dMaheswar et al.L l2010h . In the case of cores that have 
small angular sizes, we included fields containing additional 
cores that are located spatially closer and show similar radial 
velocities (obtained, for example, from CO observations), in or- 
der to have sufficient number of stars to infer their distances. 
Here we assume that the cores that are spatially closer and have 
similar velocities are located at almost the same distances. Table 
|2] summarizes the area of the fields (in degree) chosen towards 
the cores studied here, the galactic coordinates of the center of 
the fields, the total number of stars selected after applying all the 
selection criteria (cr < 0.035, SNR > 10, & (J - K,)< 0.75) and 
the number of stars classified as dwarfs by our method towards 
each field. Below we present results and discussion on individual 
clouds. In the Ay vs. d plots presented for the cores studied, the 
dash-dotted curve shows the increase in the extinction towards 
the clouds' galactic latitude a s a function of distance p roduced 
using the expressions given bv lBahcall & Soneiral(IT980h (BS80, 
hereafter). The dashed vertical line(s) is drawn at the cloud dis- 




distance (pc) 



Fig.3. Same as in Fig.|2]but for BHR 111, L328, L673, L1014 
andL1148. 



tance(s) inferred from the procedure described in ^(see Figs.|2] 
and [3]). 

3.1.1. IRAM04191 

We show the 3° x2° extinction mapFI produced bv'Dobas hi et al.l 
(2005) of the region containing IRAM04191 in Fig.E] The con- 
tours are drawn at 0.4, 0.6, 0.8 and 1.0 magnitude levels. We 
selected two fields, F 1 and F2, each covering an area of 1° x 1° 
as shown in Fig.|4] The stars, classified as dwarfs, that are used 
for determining the distance to the cloud are identified using cir- 
cles. 

In Fig. |5] we present the Ay vs. d plot for the stars from 
both Fl (filled circles) and F2 (open circles) combined together. 
The dashed vertical line is drawn at 127 pc. The jump in the 



The extinction map, covering the entire region in the galactic lati- 
tude range 1^1 < 40° derived using the optical database "Digitized Sky 
Survey I" and the traditional star-count technique, was produced in two 
angular resolutions of 6' and 18' (Dobashi et al. 2005). In this paper, we 
used the maps with 6' angular resolution. 
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Table 2. The details of the fields selected towards each cloud. 



# Area 


/ 


b 


Total 


Stars classified 


selected 






stars 


as Dwarfs 


(°) 


(°) 


(°) 






IRAM04191 










1 1.5 


178.7968 


-23.9180 


403 


277 


2 1.0 


180.2393 


-23.9180 


438 


270 






Total 


841 


547 


L1521F 










1 1.0 


172.5343 


-13.5645 


199 


155 


L l.U 


1 T 1 CI /in 

1 / i.jloy 


— 1 j.jo4j 


AAA 

444 


Jjl 


3 1.0 


171.5170 


-14.5545 


160 


120 


4 1.0 


172.5388 


-14.5545 


144 


97 


5 1.0 


171.8592 


-15.5445 


368 


243 


o l.U 


1 /U.ojz / 


— Ij.j44j 


ZSJ 


1 07 

ly / 


'7 in 
/ 1 .u 


1 7n S'^QA 
1 /U.oZy+ 


— 10. J J4J 


J JO 


JOO 


o in 

a 1.0 


1 T 1 1 n 

1 / l.oolU 


— 16.JJ43 


A10 

4j8 


oon 
ZsU 






Total 


2596 


1809 


DUD 111 

rSHKl 1 1 










i U.J 


jZ / . j4 / o 


+ i .yuuu 


^77 

JZ / 


071 
Z / i 


1 n c 
I U.5 


JZ /. 1 JZl 


+ 1 .4UUU 


Z56 


OCA 

zjU 


n c 

1 U.D 


i^n T 1 7Q 
JZ /.Zl /j 


1 n nnnn 
+U.9UUU 


7/;c 
zOj 


ZzU 






Total 


878 


741 












1 n ^ 

1 U.J 




— U.o4ol 


/4 


OZ 


2 0.5 


12.5573 


-0.5035 


328 


285 






Total 


402 


347 












1 U.J 


Afi 1 /IT/I 


— l.jljZ 


J4 


00 

zz 


2 0.5 


46.4786 


-0.8260 


108 


DO 

so 


3 0.5 


45.9856 


—0.8260 


111 


on 

89 






Total 


253 


199 


L1014 










1 n /I 
1 U.4 


m ^^77 

yz. J J / / 


— U.zo/z 


OCT 

ZJJ 


ono 

zuy 


2 0.4 


92.7318 


-0.0560 


229 


191 






Total 


482 


400 


LI 148 










1 0.7 


102.5346 


+ 15.2936 


191 


144 


2 0.7 


101.7809 


+ 15.1666 


205 


151 


3 0.7 


102.5360 


+ 16.0555 


196 


145 


4 0.7 


101.7795 


+ 15.9285 


208 


149 






Total 


800 


589 



extinction values significantly above of that expected from the 
expression of BS80 is found to occur, in both fields, at ~ 127 pc. 
Beyond this distance, the stars with high extinction are found to 
be distributed almost continuously in distance as one would ex- 
pect since the cloud will act as a dense sheet of dust layer dim- 
ming the stars shining through it. Using 541 stars that are classi- 
fied as main sequence stars, we estimated a distance of 127 + 25 
pc to the cloud IRAM04191. 

Very accurate distance determinations (< 1%) are available 
for a number of naked T Tauri stars associated with the Taurus 
molecular cloud complex from the trigonometric parallax mea- 
surements made using the Very Long Baseline Array (VLBA) 
multi-epoch observations carried out by detecting their non- 
thermal 3.6 cm radio continuum emission. These stars include, 
HDE 283572, Hubble 4, T Tau and HP Tau /G2. Th ey are found 
to be located a t distances of 128.5 + 0.6 pc (iTorres e t al., 2007), 
132j ± 0.5 pc (iTorres et a il2007h. 147 6 + 6pc Coinard et al.. 
I2007h and 161.2 ± 0.9 pc dTorres et afl.l2009h respectively. As of 
now, from the distances determined for these four stars, it is ap- 
parent that the clouds in the Taurus molecular cloud complex 
are distributed in space in distances ranging from ~ 128 pc to 
~ 161 pc. From the mean parallax obtained from the above four 
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Fig. 4. The 2° x 3° extinction map produced bv iDobashi et al.l 
(l2005l) of the region containing IRAM04191. The contours are 
drawn at 0.4, 0.6, 0.8 and 1.0 magnitude levels. The fields used 
for selecting star to determine distance are identified and la- 
belled. The stars classified as dwarfs and used for determining 
distance to the cloud are identified using open circles. 
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Fig. 5. The Ay vs. d plot for all the stars classified as dwarfs 
from the fields Fl-2 combined together towards IRAM04191. 
The dashed vertical line is drawn at 127 pc inferred from the 
procedure described in the ^(see Figl3]l. The dash-dotted curve 
represents the increase in the extinction towards the Galactic lat- 
itude of b = -24° as a function of distance produced from the 
expressions given by BS80. The error bars are shown on a few 
stars for better clarity of the points. 



stars ITorres et"al ] (12009! ) determined a mean distance of 141.2 
pc to the Taurus cloud complex. At a distance of 127 + 25 pc, 
IRAM04191 could possibly located at the front end of the Taurus 
cloud complex. In all the previous studies of IRAM04191, this 
cloud was considered to be associated with the Taurus clou d 
complex and assigned a distance of 140 pc dKenvon et al.Lll994l) . 
The distance estimate presented here is the first independent 
confirmation of the association of IRAM04191 with the Taurus 
cloud complex. 



6 



Maheswar et al.: Distances to dense cores with VeLLOs 




-90 



Arc Minutes 
Center: Longitude 171.50 Latitude 



-15.00 



Fig. 6. The ~ 3° x 4° extinction map produced by Dobas hi et alj 
(I2005h of the region containing L1521F. The contours are drawn 
at 0.5, 1.0, 1.5, 2.0, 2.5 and 3.0 magnitude levels. The fields 
used for selecting star to determine distance are identified and 
labelled. LI 52 IF is located in F3. 



3.1.2. LI 521 F 

Figure |6] shows the ~ 3° x 4° extinction map of the region con- 
taining L1521F. The contours are drawn at 0.5, 1.0, 1.5, 2.0, 
2.5 and 3.0 magnitude levels. We divided the region contain- 
ing L1521Finto eight fields, Fl-8, as shown in Fig.|6] Each field 
covers an area of 1° x 1°. The stars, classified as dwarfs, that are 
used for determining distance to the cloud are identified using 
circles. The cloud L1521F is located in the field F3. 

In Fig. 121 we present Ay vs. d plot for the stars from all the 
fields (Fl-8) combined together The dashed vertical lines are 
drawn at 1 16 and 156 pc. Though, in Fig.|7] only one component 
of dust layer is apparent, we could notice the presence of two 
dust layers in Fig.|2]as it shows a sharp rise in the mean values of 
extinction occurring at two distances, one at 1 16 pc and another 
at 156 pc. In the Ay vs. d plots of individual fields, we found 
that the presence of a dust layer at 1 16 pc is more conspicuous 
towards the fields Fl-4, than towards the fields F5-8. In Fig. [8] 
we present the Ay vs. d plot for the stars from the fields Fl- 
4 (upper panel) and F5-8 (lower panel) separately to show the 
dominance of the two components of dust layers in these two 
regions. In Fig. |2] we show the distances determined using the 
stars from Fl-4 and F5-8 independently. The two dust layers are 
found to be at distances 1 18 ± 22 and 155 ± 29 pc and 1 14 ± 21 
and 156 + 29 pc towards the fields Fl-4 and F5-8 respectively. 
The distances estimated for the two dust components using stars 
from separate sets of fields are found to be in good agreement. 
These results ascertain the fact that the technique employed here 
in estimating distances to dust layers in a direction is capable of 
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Fig. 7. The Ay vs. plot for all the stars classified as dwarfs from 
the fields Fl-8 combined together towards L1521F. The dashed 
vertical lines are drawn at distances 116 pc and 156 pc inferred 
from the procedure described in the ^ (see FigO. The dash- 
dotted curve represents the increase in the extinction towards 
the Galactic latitude of b = -14.9° as a function of distance 
produced from the expressions given by BS80. The error bars 
are shown on a few stars for better clarity of the points. 



picking them consistently even though the fields (1° x 1° each) 
used for selecting the stars are located relatively far (3-4° in this 
case) apart. 

The dust layer at ~ 156 pc seems to be present in all the 
eight fields. Also, the Ay vs. d plot produced for the sources 
from F5-8 do contain a number of stars that show high extinc- 
tion at distances smaller than ~ 156 pc. This could be due to 
the contribution from the dust component at ~ 1 14 pc. One of 
the limitations of the method is that in case of multiple jumps in 
extinction in an Ay vs. d plot, it is difficult to associate a partic- 
ular jump with any of the dust layer. In fact, this is true for all 
methods that utilize the rise in the extinction as a way to estimate 
distances. 

The above findings encouraged us to estimate distances to 
the entire Taurus molecular cloud complex. We found the pres- 
ence of four dust components at distances 85 ± 16, 114 ± 20, 
138 ± 25 and 168 ± 31 pc towards it. The results will be pre- 
sented in a forthcoming paper. The presence of the dust com- 
ponent at ~ 85 pc, though, is found to be less evident towards 
the entire cloud complex. Based on the fact that the dust layer at 
~ 1 18 pc is found to be dominant towards F3, the field in which 
L1521Fis located, we associate the jump at ~ 1 18 + 22 pc to the 
cloud L1521F. But adopting a mean value of 136 + 36 pc (error 
in this case is V2T2T292) to LI 52 IF would be more appropri- 
ate as it is closer to the distances of two sources, HDE 283572 
and Hubble 4, found to be located nearer (< 2°) to L1521F and 
are at distances 128 .5 + 0.6 pc and 132.8 ± 0.5 pc respectively 
(iTorresetaLL 120091) . 

3.1.3. BHR 111 

In Fig. |9] we present the ~ 1.5° x 1.5° extinction map of the 
region containing BHR 111 (DC 3272-1-18). The contours are 
drawn at 2, 2.1, 2.2, 2.3 an d 2.4 magnitude leve ls. The size of 
the cloud is only 10' x 4' dBourke et all Il995al) . We included 
two additional fields F2 and F3 apart from the field Fl, which 
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Fig. 8. The Ay vs. d plot for the stars classified as dwarfs from 
the fields Fl-4 in the upper panel and F5-8 in the lower panel 
towards L1521F. The dashed vertical lines are drawn at 116 pc 
and 156 pc. The dash-dotted curve represents the increase in the 
extinction towards the Galactic latitude of b = -14.9° as a func- 
tion of distance produced from the expressions given by BS80. 
The error bars are shown on a few stars for better clarity of the 
points. 



contains BHR 111, as shown in Fig.|9] Unfortunately, no velocity 
information is available on the regions selected except for BHR 
111 dBourke et al.L il995al) . The stars, classified as dwarfs, that 
are used for determining the cloud distance are identified using 
circles. 

In Fig. [TO]we present the Ay vs. d plot for the stars from all 
the fields (Fl-3) combined together. The dashed vertical line is 
drawn at 355 pc. There are a number of stars that exhibit rel- 
atively high extinction values at distances smaller than 355 pc. 
Inspection of Fig.[TT] which presents the Ay vs. d plot for stars 
from the individual fields separately, shows that the jump in the 
extinction values consistently occurs at or beyond 355 pc. The 
presence of a foreground dust layer is suspected. On the basis 
of 741 stars classified as dwarfs by the method, we estimate a 
distance of 355 ± 65 pc to BHR 111. 

Distance to BHR 111 was first determined by iBourke et alJ 

(Il99 5b) from a plot of stellar reddening of stars, selected from 
a circular region about the object, as against their distances cor- 
rected for reddening. They selected stars that had MK spectral 
types and colours available. They estimated a distance of 250 pc 
to BHR 111. Using the same method. lRacca et al. ( 2009) also es- 
timated a distance of 250 pc to BHR 1 1 l. fBourke et aL d 1995b) 
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Fig. 9. The ~ 1.5 ° x 1.5° extinction map produced by 
Do bashi et al.1 (l2005h of the region containing BHR 111. The 
contours are drawn at 2, 2.1, 2.2, 2.3 and 2.4 magnitude lev- 
els. The fields used for selecting star to determine distance are 
identified and labelled. The stars classified as dwarfs and used 
for determining distance to the cloud are identified using open 
circles. 
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Fig. 10. The Ay vs. d plot for all the stars classified as dwarfs 
from the fields Fl-3 combined together towards BHR 111. The 
dashed vertical line is drawn at 355 pc inferred from the pro- 
cedure described in the ^ (see FigO. The dash-dotted curve 
represents the increase in the extinction towards the Galactic lat- 
itude of b - 1 .824° as a function of distance produced from the 
expressions given by BS80. The error bars are shown on a few 
stars for better clarity of the points. 
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Fig. 11. The Ay vs. c/ plot for the stars classified as dwarfs from 
the individual fields Fl, F2, & F3 towards BHR 111. The dashed 
vertical line is drawn at 355 pc inferred from the procedure de- 
scribed in the ^(see FigO. The dash-dotted curve represents 
the increase in the extinction towards the Galactic latitude of 
b = 1.824° as a function of distance produced from the expres- 
sions given by BS80. The error bars are shown on a few stars for 
better clarity of the points. 



used a search radius of 5° about the cloud to obtain stars with 
MK spectral types and colours. If not enough stars were avail- 
able within 5°, the y increased the search radius to 7.5° or 10°. 
iRacca et al.1 (l2009l) ' used a search radius of 3° about the cloud to 
select the stars. In Fig.|9l it can be noticed that the cloud is not 
isolated but is surrounded by diffused dust clouds. In Fig.fTTIwe 
noticed the evidence of more intervening dust layers towards the 
line-of-sight which is more conspicuous towards F 3. It could be 
poss ible that the large search radius used by both Bourke et alj 
d 1995 b) and Racca et al. (2009) might have picked the nearer 
dust layer possibly located at 250 pc. 



3.1.4. L328 

In Fig. [12] we present the 1° x 1° extinction map of the region 
containing L328. The contours are drawn at 1.5, 2.0, 2.5, 3.0, 
3.5 and 4.0 magnitude levels. We divided the region containing 
L328 into two fields, Fl and F2, as shown in Fig. [12] The angular 
extend of L328 (~ 1') is found to be too small to divide the field 
containing the cloud further into smaller sub-fields with suffi- 
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Fig. 12 . The ~ 1 ° x 1 ° extinction map produced bv lDobashi et al.l 
(I2OO5I) of the region containing L328. The contours are drawn 
at 1.5, 2.0, 2.5, 3.0, 3.5 and 4.0 magnitude levels. The fields 
used for selecting star to determine distance are identified and 
labelled. The stars classified as dwarfs using the near-IR pho- 
tometry in each field are shown using filled circles. The fields 
Fl and F2 cover an area of 0.3° x 0.3° and 0.5° x 0.5° respec- 
tively. 
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Fig. 13. The Ay vs. d plot for all the stars classified as dwarfs 
from the fields Fl (open circles) & F2 (filled circles) combined 
together towards L328. The dashed vertical line is drawn at 217 
pc inferred from the procedure described in the ^(see Fig[3]l. 
The dash-dotted curve represents the increase in the extinction 
towards the Galactic latitude ofb- -0.829° as a function of dis- 
tance produced from the expressions given bv lBahcall & Soneiral 
(I198Q) . The error bars are shown on a few stars for better clarity 
of the points. 



cient number of stars for the analysis. Therefore we included the 
field containing an another molecular core, L323, located in the 
close vicinity (~ 20') of L328. Both the cores, L328 and L323, 
are found to have similar radial velocities of 6.6 and 6.4 km s ' 
respectively dClemens &BarvainisL [19881) . The locations of the 
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Fig. 14. The 1.5° x 1.5° extinction map produced by 
iDobashi et alj ( l2005h of the region containing L673-7. The con- 
tours are drawn at 4.0 and 5.0 magnitude levels. The fields used 
for selecting star to determine distance are identified and la- 
belled. The stars classified as dwarfs using the near-lR photome- 
try in each field are shown using filled circles. Each field covered 
an area is 0.5° x 0.5°. 



cores and the fields chosen are identified and labelled in the Fig. 
[T2l The stars, classified as dwarfs, that are used to estimate the 
core distance are identified using circles. 

In Fig. [13] we show the Ay^s,- d plot for the stars from all the 
fields (Fl-2) combined together. The stars from Fl and F2 are 
represented using open and filled circles respectively. The field 
containing L328 is small when compared to F2 due to the small 
angular extend of the cloud. The stars showing high extinction 
at/or beyond 217 pc are present in both the fields. There exist 
two stars with Ay > I magnitude at distances smaller than 217 
pc. Both stars are found to be located towards F2 at different dis- 
tances. On the basis of 347 sources classified as main sequence 
stars from the two fields, we estimated a distance of 217 ± 30 pc 
to bothL328 and L323. 

Earlier estimates of distances to L323 (and L328) were by 

iBok & McCarthv (1974) who assumed a distance of 200 pc 
based on the number of stars brighter than nipg = 21 mag pro- 
jected on the cloud in their photographic plate. Later, several 
authors have assigned 200 pc assuming that the clouds, L328 
and L323, are located in between the Ophiuchus cloud and the 
Aquila Rift (e.g., [Le e et al., 2009). Using photometry in Vilnius 
system. [Straizvs etalj (12003 ) estimated a distance of 225 ± 55 pc 
to the front edge of the Aquila Rift with a possible depth of 80 
pc to the cloud complex. If L328 and L323 are physically asso- 
ciated with the Aquila Rift, then possibly, the clouds are located 
to the front edge of it. 

3.1.5. L673-7 

In Fig. [14] we present the 1.5° X 1.5° extinction map towards 
the direction of L673-7. The contours are drawn at 4.0 and 5.0 
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Fig. 15. The Ay vs. d plot for all the stars classified as dwarfs 
from the fields Fl, F2 & F3 combined together towards L673. 
The dashed vertical line is drawn at 240 pc inferred from the 
procedure described in the ^(see Fig[3]l. The dash-dotted curve 
represents the increase in the extinction towards the Galactic lat- 
itude of ^7 = -1.33° as a function of distance produced from the 
expressions given by BS80. The error bars are shown on a few 
stars for better clarity of the points. 



magnitude levels. We divided the region containing L673-7 into 
three fields, Fl, F2 and F3, as shown in Fig. [14] L673-7 is 
part of the L6 73 cloud complex from the Lvnds (1962) catalog. 
iLee & MversI tl999 ) identified 1 1 cores within this complex in- 
cluding L673-7. The radial velocity of L6 73 complex in CO line 
is shown to range from6.71 - 7.30 km s ' (iKislvakov & GordonL 
1983)- Some more clouds are found to be associated with the re- 
gion s elected around L673 namel y, L675, L676 from th e lLvndsl 
(11962!) catalog and CB188 from Iciemens & Barvairiisl (Il988h . 
These three clouds also sho w similar radial velocities of 7.5, 7.6 
and 7. 1 km s ' respectively (Clemens & B aryainislll988h as that 
of L673 complex. We therefore assume that the regions selected 
around L673 are all associated and included in our analysis. The 
stars, classified as dwarfs, that are used for determining the cloud 
distance are identified using circles. 

In Fig. [15] we present the Ay vs. d plot for the stars from 
all the fields (Fl-3) combined together A jump in the extinction 
values significantly above the values expected from the expres- 
sion of BS80 is evident at/or beyond ~ 240 pc. The stars showing 
high extinction are distributed almost uniformly beyond ~ 240 
pc. In Fig. [16] we show the Ay vs. d plot for the stars from the 
individual fields Fl, F2 and F3. In all the three fields, the jump 
in the extinction is evidently occurring at ~ 240 pc. Using 199 
sources classified as main sequence stars from the three fields, 
we estimated a distance of 240 + 45 pc to L673-7. 

Previo us determination of dist ances to L673-7 were highly 
uncertain. iHerbig & Jones! (!l983l) used a distance of 300 pc, 
based on proper motion studies. Others assumed a distance of 
either 200 pc or 300 pc in th ei r studies (e. g., Hilton & LahullaL 
!1995L !Dunham et all !2010al) . !Edwards & Snell (,1982h used a 
distance of 150 pc based on the assumption that L673 is a part 
of Gould's belt and characterized it as the minimum distance to 
the cloud. Because L673, L675, L676 and CB188 show similar 
radial velocities, we assign a distance of 240 + 45 pc to all these 
clouds. 
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Fig. 17. The 2° x 2° extinction map produced by (iDobashi et al.Ll2005l) of the region containing LDN 1 148. The contours are drawn 
at 0.7, 1.0 and 2.0 magnitude levels. The fields used for selecting star to determine distance are identified and labelled. Each field is 
of0.7°x0.7°. 



3.1.6. L1148 

The fields used to determine the distance to L1148 are shown 
on the extinction map in Fig. [17] The contours are drawn at 0.7, 
1.0 and 2.0 magnitude levels. We divided the region containing 
LI 147/48 into four fields, Fl, F2, F3 and F4, as identified on Fig. 
[TtI The stars, classified as dwarfs, that were used to determine 
the cloud distance are identified using circles. 

In Fig. [T8]we present the Ay vs. d plot for the stars from all 
the fields (Fl-4) combined together Ajumpin the extinction val- 
ues significantly above the values expected from the expression 
of BS80 is evident at 168 pc. But if we compare the Ay vs. c/ plots 
of 1RAM04191 and LI 52 IF, which are both located within 200 
pc distance, with that of LI 148, we can notice that the stars with 
high extinction towards LI 148 are not distributed continuously 
beyond 168 pc. The angular size of the cloud (inferred from the 
extinction map shown in Fig.fTTTi is large enough to have shown 
a continuous distribution of stars with high extinction beyond 
168 pc if the cloud would have been located at this distance. The 
lack of such an effect is clearly seen in Fig. [T9l where we show 
the Ay vs. d plot for the stars from the individual fields Fl (open 
circles) and F2 (filled cicles) in the upper panel and F3 (filled 
circles) and F4 (open circles) in the lower panel. Only few stars 
are seen with high extinction especially towards the fields F3 and 
F4. On the contrary, a wall of high extinction stars are seen be- 
yond ~ 300 pc towards all the four fields. Using the stars from 
F3 and F4 alone we estimated a distance of 301 + 55 pc to this 
dust component. 

The LI 148 core is located in a cloud complex towards 
the direction of Cepheus. The most widely quoted distance to 
LI 147/48/55 is 325 + 13 pc (Straizys et al., 1992). They used the 
Vilnius photometric system to classify the stars and also to get 
the interstellar reddening. They obtained a distance of 325 + 13 



pc to LI 147/1 158 by taking an average of ten stars sh owing 
Ay > .45 which were distributed in the range 240-380 pc. lSnelll 
(fl98Th obtained a distance of 250 pc to LI 147 based on the Ay 
vs. d plot produced using the stars from a region of 20° x 20° 
centred on the cloud t hat h ad both optical photometry and MK 
spectral classification. lKurJ (11998) showed the presence of three 
dust components at characteristic distances: 200, 300, and 450 
pc, based on cumulative distribution of field star distance mod- 
uli in Wolf diagrams. A complex velocity structure was reported 
bv Hariunpaa et al.. (.1991) towards the direction of LI 147/48, 
LI 152, LI 157 and LI 155/58 cloud complex. They found the 
presence of two distinct velocity components in '^CO and C'^O 
line profiles. The blueshifted component {Vlsr ~ 1.5 km s"') 
was found to be limited to smaller extend whereas the redshifted 
component {Vlsr ~ 2 - 4 km s"') was found to be extended 
over the whole cloud area. Based on the velocity structure found 
over the whole cloud, they suggested the presence of two cloud 
components along the line of sight. 

The radial velocities measured for the cores LI 148, LI 155C- 
1 and L1155C-2 of the cloud comp lex are fo und to be 2.56, 
2.58 and 1.35 km s"' respectively (Lee et al.L ,2004). Another 
cloud, LI 167/1 174, which al s o showed similar radi a l velo city 
(2.67 km s-\ ICaselli et all llool lYonekura et all [T997h as 
that of LI 147/48 is found to be located at ~ 2° south-east of 
LI 147/48/55 cloud complex (see Fig. [TTll. IStraizvs et all (Il992h 
estimated a distance of 288 ± 25 pc to the LI 167/1 174 again by 
taking an average distance of four considerably reddened stars. 
They pr eferred a distanc e of 275 pc to HD 200775, a Herbig 
Be star dThe et all [19941) found associated with NGC 7023 and 
responsible for the reflection nebulosity, by assuming it to be 
a B3Ve star Maheswar et al. (2010), using the same method as 
employed in this work, have estimated a distance of 408 ± 76 
pc to NGC 7023 (consistent with the Hipparcos parallax dis- 
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Fig. 16. The Ay vs. d plot for the stars classified as dwarfs 
from the individual fields FI (upper), F2 (middle) and F3 (lower 
panel) towards L673. The dashed vertical line is drawn at 240 
pc inferred from the procedure described in the ^(see FigO. 
The dash-dotted curve represents the increase in the extinction 
towards the Galactic latitude ofb- -1.33° as a function of dis- 
tance produced from the expressions given by BS80. The error 
bars are shown on a few stars for better clarity of the points. 
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Fig. 18. The Ay vs. d plot for all the stars classified as dwarfs 
from the fields FI, F2, F3, & F4 combined together towards 
L1148. The dashed vertical line is drawn at 168 pc inferred 
from the procedure described in the ^ (see FigO. The dash- 
dotted curve represents the increase in the extinction towards the 
Galactic latitude of b - -1-15.3° as a function of distance pro- 
duced from the expressions given by BS80. The error bars are 
not shown on all the stars for better clarity. The sources shown 
using filled circles in red are those identified with black squares 
in Fig. [13 



tance to this cloud by selecting stars from a region of 1° x 1° (see 
Fig. 117b centred on the cloud. In Fig.|2Tl we show the Ay vs. d 
plot for the stars towards L1122. A sudden jump in the extinc- 
tion is visible at 220 pc. We assigned the distance of the first star 
that showed Ay > 1 as the distance to this cloud. It should be no- 
ticed that the dust components which were seen at 310 and 408 
pc towards LI 167/1 174 are absent towards LI 122. The location 
of this cloud at a smaller distance proves that there exist dust 
components that are as close as 200 or even less as we found to- 
wards LI 148. We therefore believe that the dust component seen 
towards LI 147/48 cloud complex is a foreground layer and that 
the actual layer is associated with the dust component inferred 
at 301 pc. 



tance of429|'^Sc to HD 200775. lvan den Ancker etaTlllQgSt 
iBertout et all 119991) . In the direction of NGC 7023, they also 
found evidence for two additional dust components one at ~ 200 
pc and another at 305 pc w hich were in good agreement with 
those inferred bv iKunI (11998) . In Fig. |20]we show the vs. 
d plot for the stars from the fields F5 and F6 which contain 
LI 167/1 174 complex (see Fig.fTTll. Step-like features at 310 and 
408 pc are clearly visible in the plots. The component at 408 pc 
seems to be more conspicuous towards F5 where LI 174 cloud is 
located. A continuous distribution of stars with high extinction is 
apparent beyond 408 pc. HD 200775 is probably associated with 
this cloud component. Towards F6, the dust component found at 
310 pc seems to be dominant. There are a number of stars both 
towards F5 and F6 that show high extinction even at distances 
smaller that 310 pc. Presence of a third dust layer in the line of 
sight could be the reason. 

About 2° south-west of the LI 147/48 cloud complex, there 
exist one more Lynds cloud , LI 122. This c l oud s howed a ra- 
dial velocity of 4.8 km s ' (lYonekura et al.L Il997h . No earlier 
distance estimates are available for the cloud. We estimated dis- 



3.1.7. L1014 

In Fig. |22] we identify the fields used for selecting stars on the 
1.2° X 1.2° extinction map of L1014. The contours are drawn at 
1.0, 1.5 and 2.0 magnitude levels. The stars, classified as dwarfs, 
that are used for determining the cloud distance are identified 
using circles. Each field covers an area of 0.4° x0.4°. The angular 
extend of L1014-2, which contains the VeLLO, L1014-IRS, is 
only a: 2'. Therefore we included another cloud L1014-1 which 
is located ~ 10' north-east of L1014-2. The cloud L1014-1 is 
better known as B362. Both L1014-1 and L1014-2 are found to 
have similar radial velocities (4.0 and 4.2 km i"' respectively, 
Lee et al., 1999; Crapsi et al., 2005). 

In Fig. |23] we present the Ay vs. d plot for the stars from all 
the fields FI (filled circles) and F2 (open circles) combined to- 
gether In Fig. [3] the steep rise in the mean values of extinction 
are found to happen at distances 221 and 258 pc. But from Fig. 
l23] it is evident that the steep rise in the extinction at 221 pc is 
occurring mainly due to two stars that are located in the field, 
F2 alone. On the contrary, the steep rise in the mean value of 
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Fig. 19. The Ay vs. d plot for the stars classified as dwarfs from 
the individual fields Fl (open circles) and F2 (filled circles) in 
upper panel and F3 (filled circles) and F4 (open circles) in lower 
panel towards LI 148. The dashed vertical lines are drawn at 168 
pc and 301 pc inferred from the procedure described in the ^ 
(see Figl3]l. The dash-dotted curve represents the increase in the 
extinction towards the Galactic latitude of b = +15.3° as a func- 
tion of distance produced from the expressions given by BS80. 
The error bars are shown on a few stars for better clarity of the 
points. 



Fig. 20. The Ay vs. d plot for the stars classified as dwarfs from 
the individual fields F5 (upper panel) and F6 (lower panel) to- 
wards NGC 7023. The dashed vertical lines are drawn at 168 
pc, 310 pc and 408 pc inferred from the procedure described in 
the ^(see Fig|3]l. The dash-dotted curve represents the increase 
in the extinction towards the Galactic latitude of b = -h14.1° as 
a function of distance produced from the expressions given by 
BS80. The error bars are shown on a few stars for better clarity 
of the points. 



extinction at 258 pc (see Fig. l23T l is caused due to the stars from 
both fields, Fl and F2 (filled and open circles). Therefore we as- 
sociate the jump in the extinction at 258 pc to the presence of the 
cores L1014-1/2. Also, beyond this distance the high extinction 
stars are distributed almost regularly. However, we found only 
two stars with low extinction values lying foreground to L1014- 
1/2 that are essential to have a complete census of the dust layers 
in the direction. In order to include more cores in the vicinity 
of L1014-1/2 that share similar radial velocities, we looked at 
the lar ge-scale maps of the clouds produced by i Dobashi et alJ 
(11994 based on '"'CO (with a 2.7' angular resolution) survey 
towards the Cygnus region. L 1014-1, the cloud which was de- 
tected in their survey, is found to be isolated with no cloud com- 
ponents, at least, within a radius of 2° in the velocity range of 
3 < Vlsr < 4.5 km i"' (see Fig. 5/ of their publication). 

iDame & Thaddeusl(ll985h have assigned a dis tance of 800 pc 
to L1014 assuming it to be part of Cyg OB7 (see lDobashi et al.L 
1994). Ho et al. ( 1978) adopted a distance of 200 pc for B362. 
Morita et all (|2006|) assigned a distance of 400-900 pc to L1014 
by assuming possible association of identified T Tauri stars with 



the cloud. No reliable estimate of distance to L 1014-2 is avail- 
able in the literature. In most of the studies carried out on L1014- 
2, the author s have adopted a distance of 200 pc to it (e.g., 
lLee&Mverill999h bv assuming that both B362 and L1014-2 
are located at same distance. Using a total of 400 sources classi- 
fied as main sequence stars from the two fields, we estimated a 
distance of 258 ± 50 pc to L1014. 

3.2. Are VeLLOs really VeLLOs? 

Using the distances estimated to the seven cores in this work we 
re-evaluate the internal luminosities of the VeLLO candidates. In 
Table|3]we list the distances used in earlier studies and the inter- 
nal luminosities calculated using those distances of the VeLLOs 
in columns 2 and 3 respectively. The re-evaluated internal lumi- 
nosities of the VeLLOs are given in column 6 of the Table [3] 
The previously determined internal luminosities of the VeLLOs 
in the cores L328, L673-7 and L1014 were calculated using dis- 
tances that were highly uncertain. However, we find that these 
distances are very close to the values we have determined in this 
work. As we can see from the Table |3] that the re-estimation of 
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Fig. 21. The Ay vs. d plot for all the stars classified as dwarfs 
towards LI 122. The dashed vertical lines are drawn at 220, 301 
and 408 pc. The dash-dotted curve represents the increase in the 
extinction towards the Galactic latitude of b = +14.8° as a func- 
tion of distance produced from the expressions given by BS80. 
The error bars are shown on a few stars for better clarity of the 
points. 




Arc Minutes 
Center: Longitude 92.47 Latitude -0.20 

Fig. 22. The 1.2° x 1.2° extinction map produced by 
[Dobashi et al.l (1200 5') of the region containing L1014. The con- 
tours are drawn at 1.0, 1.5 and 2.0 magnitude levels. The fields 
used for selecting star to determine distance are identified and 
labelled. The stars classified as dwarfs using the near-IR pho- 
tometry in each field are shown using filled circles. Each field is 
of 0.4° X 0.4°. 



internal luminosities of six VeLLO candidates with our reliable 
distance measurements confirm them to be consistent with the 
definition of the VeLLOs. The VeLLO associated with L1014- 
2 core is found to be L,„, = 0.15 Lq. However, this source is 
still interesting and could be called as VeLLO-like since the lu- 
minosity is still an order of magnitude less than the accretion 



Fig. 23. The Ay vs. d plot for all the stars classified as dwarfs 
from the fields Fl & F2 combined together towards L1014. 
The dashed vertical lines are drawn at 221 and 258 pc inferred 
from the procedure described in the ^ (see Figl3]l. The dash- 
dotted curve represents the increase in the extinction towards the 
Galactic latitude of b = -0.251° as a function of distance pro- 
duced from the expressions given by BS80. The error bars are 
shown on a few stars for better clarity of the points. 

Table 3. On the status of the VeLLOs. 



Cores d' Lj,,, Ref. d"^ L, 



(pc) (Le) (pc) (Le) 



IRAM04191 


140 


0.08 


1 


127 ± 


25 


0.07 


L1521F 


140 


0.05 


2 


136 ± 


36 


0.05 


BHR 111 


250 


0.04 


3 


355 + 


65 


0.08 


L328 


200 


0.06 


4 


217 + 


30 


0.07 


L673-7 


300 


0.04 


5 


240 ± 


45 


0.03 


L1014 


200 


0.09 


6 


258 ± 


50 


0.15 


L1148 


325 


0.10 


7 


301 ± 


55 


0.09 



' Distances used by authors to evaluate internal luminosities of the 

VeLLOs p rior to t h is wor k. ^Distance estimates from this work. 

1. Dunha m etalJ (l2006h: 2. [Bourke et al. (2006); 3. Dunh am et al. 
(2008); 4. Lee et al.' (200^); 5. iPunham et al.i l i2010ai) ; 6. lYoung et al. 
(2004); 7. Kauffmann et al. (20051) ; 

luminosity of Lace ~ 1-7 estimated for a 0.08 Mq protostar 
using an accretion rate of ~ 2 x 10'^ yr" ' (the rate predicted 
by the standard model; I Shu et all 1 19871) and 3 Rq stellar radius 
in the expression. Lace = 3.13 x 10^ MMgcdR Lq. 

4. Conclusions 

We estimated distances to seven dark clouds IRAM04191, 
L1521F BHRlll, L328, L673, L1014, and LI 148 which are 
found to harbour VeLLO candidates, discovered in the Spitzer 
Legacy program — From Molecular Cores to Planet Forming 
Disks, using near-IR photometric method. In this method, first 
the 2MASS JHR, photometry are used to produce the (J - H) 
and (H - K,) colours of the stars projected on the field con- 
taining each cloud. These observed colours are then dereddened 
simultaneously using trial values of Ay and a normal interstel- 
lar extinction law. The best fit of the dereddened colours to the 
intrinsic colours giving a minimum value of then yields the 
corresponding spectral type and Ay for the star. The main se- 
quence stars, thus classified, are then utilized in an Ay versus 
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distance plot to bracket the cloud distances. The typical error in 
the estimation of distances to the clouds are found to be ~ 18%. 
Our distance estimates are: 127 ± 25 pc (IRAM04191), 136 ± 36 
pc (L1521F), 355±65 pc (BHRl 1 1), 217±30pc (L328), 240±45 
pc (L673), 258 + 50 pc (L1014), and 301+55 pc (LI 148). Using 
these distance estimates, we re-evaluated the internal luminosi- 
ties of the VeLLO candidates discovered in these cores. Except 
L1014-IRS (Li„, = 0.15 Lo), all other VeLLO candidates are 
found to be consistent with the definition of a VeLLO. 
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